Quantum Entanglement and Detection of Topological Order in Numerics 
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"Topological ordered" phases such as gapped quantum spin-liquids and fractional quantum Hall 
states possess ground state degeneracy on a torus. We show that the topological nature of this 
degeneracy has interesting consequences for the entanglement structure of the degenerate ground 
states. This leads to a simple method for detecting topological order, which is ready-made for 
numerical schemes such as density matrix renormalization group. The method is valid for phases 
that do or do not possess edge states alike. We demonstrate it by calculating the entanglement 
spectrum and the topological entanglement entropy for a chiral spin-liquid state on a 18 site system. 

PACS numbers: 



Entanglement properties of the ground state wavefunc- 
tions have lead to a significant progress in our numeri- 
cal ability to tackle strongly correlated systems and have 
also shed light on the conceptual foundations of such sys- 
tems In the realm of gapped quantum phases, one 
of the most interesting discoveries has been the fact that 
the so-called "topologically ordered phases" that 
is, phases that have anyonic excitations in the bulk, ex- 
hibit long-range entanglement in their ground state wave- 
function One manifestation of the long range en- 
tanglement is a non-zero Topological Entanglement En- 
tropy (TEE) jtopo the entanglement entropy S 
corresponding to the ground state of a two dimensional 
topologically ordered system for a disk-shaped region A 
of linear length L with smooth boundaries scales as 



S{L) = L- jtopo + 0{l/L) 



(1) 



The number ^topo is a robust property of the gapped 
phase itself and equals \og{D) where D is the total quan- 
tum dimension [13, it does not change if the disk- 
shaped region A is deformed in any way, while maintain- 
ing its topology and keeping the boundaries smooth. This 
expression has been verified numerically for various topo- 
logically order states 15 - 2^ ■ We also note that when the 
topological ordered state has gapless edge modes, an ele- 
gant and powerful way to identify topological order is to 
study the entanglement spectrum [2l| i.e. the spectrum 
of eigenvalues of the reduced density matrix pA associ- 
ated with region A. 

As an another application of the entanglement struc- 
ture of ground state wavefunctions, the numerical 
methods such as density matrix renormalization group 
(DMRG) 0, [2^ and multiscale entanglement renor- 
malization ansatz (MERA) implement a real-space 
renormalization group that successively coarse-grain the 
system on the basis of the entanglement between a sub- 
region and the rest of the system. In particular, and 
what suffices to note for our purposes, DMRG involves 
dividing the total system, often a cylinder or a torus, 
into two halves (say A and B) that share a boundary of 
length L, and as an output, yields a Schmidt decomposed 
variationally optimized ground state wavefunction | : 



(2) 



Here {Aq,} and {|ua), \va)} are the Schmidt eigenvalues 
and eigenvectors respectively. This decomposition not 
only provides one with the correlations of various local 
operators, it also readily yields the Renyi entropies Sn — 
— ln(^^ AJ^) as well as the von Neumann entropy 
SvN ~ ^ 'Ylia In Aq associated with the bipartition. 

Juxtaposing these two developments, namely, a non- 
zero "ftopo for topologically ordered phases and the ad- 
vent of methods such as DMRG, one would think that 
entanglement based numerical techniques could be used 
to search for new topological phases, especially when 
quantum Monte Carlo can not be performed because 
of fermionic sign problem. Indeed, search for quantum 
spin-liquids using methods that exploit entanglement is 
an active area of research 25-23. In particular, Ref. 
[29I has suggested that the kagome lattice antiferromag- 
net Heisenberg model (KAH) is a gapped quantum spin- 
liquid. The authors of Ref. [l^ also found an almost 
degenerate set of two lowest lying states on a cylinder, 
which is consistent with Z2 spin liquid ground state. 
However, to establish topological order, it is desirable 
to extract a number such as jtopo, especially because one 
already has direct access to the entanglement spectrum 
{Aq} in DMRG. Unfortunately, there are at least a few 
practical obstructions for using DMRG to extract "ftopo- 
Eqn[T] suggests that one may be able to extract ^topo by 
calculating S" as a function of system size L for L ^ ^, 
the correlation length and extrapolating it to L — )■ to 
obtain jtopo- In practice, this is rather difficult to im- 
plement because the two-dimensional DMRG is typically 
limited to cylinders or strips of width L < 10 lattice spac- 
ings and the finite size effects incurred while extrapolat- 
ing could be significant. A much more accurate method 
to extract TEE is to divide the total system into regions 
A, B, C and D and calculate the following tripartite en- 
tanglement associated with regions A, Band C 



■ itopo = Sa + Sb + Sc — Sab — Sbc — Sca + Sabc (3) 
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The above expression cancels out the leading 'area-law' 
term and indeed yields an expression for ^topo that is 
accurate upto 0(1/L). However, calculating expression 
such as EqI3] within DMRG is rather challenging and has 
not been attempted so far to the best of our knowledge. 
Furthermore, since this method involves dividing the sys- 
tem into four parts (compared to two in a conventional 
DMRG setting) and calculating their entanglement, the 
finite size effects can spoil the extraction of ^topo as they 
will be determined by the smallest of the regions A-D. 

In this paper, we devise a technique that can be used 
to cancel out the non-universal leading area law term to 
extract ^topo, without resorting to the calculation of an 
expression such as Eq\^ and which can be implemented 
directly within DMRG or its variants. To begin with, we 
note that apart from the presence of a non-zero jtopo, one 
other important consequence of topological order is that 
it results in a degenerate set of ground states, on a 
non-trivial manifold such as a torus or a cylinder {i = 1- 
iV where is the ground state degeneracy). The set of 
degenerate ground states have the property that in the 
thermodynamic limit, the expectation value of any local 
operator O has the form {£,i\0\£,j) = OSij, where is O is a 
number that is independent of the index i, j. In a recent 
paper [l9l | it was shown that one can extract the braiding 
and statistics of anyonic quasiparticles from the TEE of 
the degenerate set of ground states. We now show that 
the generic properties of the entanglement spectrum of 
degenerate ground states impose strong constraints on 
the full entanglement entropy itself. This could allow 
one to detect topological order just using the bipartite 
entanglement in a finite sized system, provided the cor- 
relation length of local operators ^ ^ L, the system size, 
and the splitting of ground state degeneracy is negligible 
compared with the gap above the ground states. 

We assume that the entanglement entropy of a gapped 
phase in two dimensions may be separated as: 



S — Slocal + 



topo 



(4) 



Here Siocai is obtained by patching local contributions 
from the entangling surface and therefore admits a cur- 
vature expansion [30|. On the other hand, Stopo cannot 
be obtained in the same way, and hence is appropriately 
identified as the topological part of the entanglement en- 
tropy —jtopo- We note that Eq|4] is not merely a def- 
inition of Stopo and it will be violated if, for example, 
there are terms of the form Siocai x Stopo- To leading 
order, Eq|4] is known to hold for all gapped phases in 
two dimensions and this will be sufficient for our pur- 
poses. Since the degenerate set of ground states are lo- 
cally indistinguishable and possess the same correlation 
functions for all local operators, Siocai is same for all de- 
generate ground states. However, Stopoiogicai in general 
does depend on the ground state, if the boundary of the 
region A is non-contractible |31| . In the following, we 



will restrict ourselves to states that only have excitations 
with abelian statistics. 

Let us therefore consider a topologically ordered state 
that has N degenerate ground states on a torus or cylin- 
der. We are interested in the bipartite entanglement en- 
tropy of the ground states for an entanglement cut sim- 
ilar to the one shown in FiglT] that divides the system 
into two halves. As introduced in Ref.[3l, there exist a 
special set of ground states, dubbed 'minimum entropy 
states' (MES) jS^), i = for which the entanglement 
entropy is minimum (or in other words, the magnitude 
of the TEE is maximum). Each MES |Si) is associated 
with a distinct particle type threading perpendicular to 
the entanglement cut. As a concrete example, in the case 
of quantum Hall systems on a torus, |Si) 's are eigenstates 
of the Wilson loop operator of the Chern-Simons gauge 
field that threads flux normal to the entanglement cut; 
while for the discrete gauge theories on a torus, |Si)'s 
are the eigenstates of the operators that thread an elec- 
tric and magnetic field normal to the cut 3 

In an 

abelian theory, one is free to permute the particle labels 
arbitrarily; therefore the Renyi TEE, and hence the full 
entanglement spectrum is same for all Furthermore, 
the same considerations imply that the Schmidt states for 
different MESs arc orthogonal to each other: 



(5) 



where 



SijSaB- The factor of 



Sap follows from the definition of Schmidt decomposi- 
tion while the factor of 5ij follows from the fact that 
each Schmidt state carries the distinct quantum 

number associated with the operator that measures the 
quasiparticle type threading the entanglement cut (here, 
i). It follows from Eql5] that a generic wavefunction 
I*) = E»c»|S,) = EiaCiv^K) ® \vi) is already 
Schmidt decomposed when expressed in terms of the 
MES [s^l- Therefore, the n'th Renyi entropy correspond- 
ing to the entanglement bipartition such as that shown 
in FiglU for the wavefunction Y^) is. 



Sn — 



In(E.AS) ln(EJc.P") 



local 



bo HD) - 



n — 1 

HE. 



n - 1 



(6) 



In deriving EqnlHlwe have used the fact that Sn,iocai 
is same for all degenerate ground states and that the 
TEE for MESs equals bo ln(D) where bo is the number 
of connected components of the boundary of region A 
(for example, 6o = 1 (2) when the overall geometry is a 
cylinder (torus) js^])- Eqs. [5] and |6] can be used to es- 
tablish topological order and extract TEE from bipartite 
entanglement entropy. We note that EqlS] agrees with the 
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FIG. 1: Bipartition for the CSL state for which Schmidt de- 
composition is performed. Regions A and B each are of size 
3x3 and the spins are located at the vertices of the shown 
lattice. The wavefunction |0, 0) and |0, tt) are obtained by im- 
posing periodic-periodic and periodic-antiperiodic boundary 
conditions respectively for the slave fermions along the x, y 
directions. 




topological entanglement entropy of Clicrn-Simons gauge 



theory calculated in Ref.jSl 



We emphasize that the dependence of Sn on the ground 
state wavefunction in Eq\^results solely from the fact that 
the boundary of region A is non-contractible and there- 
fore, is of topological nature. For contractible bound- 
aries, all topologically degenerate ground states will have 
the same entanglement entropy. This, along with EqjSJ 
provides a sharp distinction between a topologically or- 
dered phase and a short-range entangled phase where the 
ground state degeneracy may result from spontaneous 
symmetry breaking (ssj . 

As a demonstration of the ideas developed above, 
we first note that the EqsH] and [6] are readily satisfied 
by fixed point wavefunctions that have zero correlation 
leng th, such as the ground states of the toric code model 
[lol| on a torus or cylinder [l^. As a non-trivial test, 
we consider the entanglement spectrum of an SU (2) spin 
symmetric Chiral Spin Liquid (CSL) state [35l - l37| on a 
rectangular lattice of size 6x3 (FiglT]). Note that, un- 
like toric code model [llj, we do not have access to the 
MES apriori. The topological order in the CSL state has 
already been established numerically in Ref.[l3| by cal- 
culating 7topo and verifying that it matches the theoret- 
ical value 7topo = log(2)/2. The ground state wavefunc- 
tion for the CSL state are obtained using a slave particle 
construction where one writes the spins as bilinears in 
fermions S ~ \ fl\(Ann' fy , and assumes a chiral d-wave 
state for the fermions 35, ^61 with a mean-field BCS 



gap A. To obtain a spin wave function, one Gutzwiller 
projects the d^^ _y2 -|- idxy superconductor to one fermion 
per site. The degenerate ground states are obtained by 
putting the system on a torus and choosing periodic or 
anti-periodic boundary conditions for the slave fermions 
along the two torus directions 0, 0, 36 1. We denote 



the ground states by the mean field fluxes in the x and ij 
directions as |<^i,y'2), 'fi,2 = 0,7r. The two fold degener- 
acy of the CSL implies that only two of the four ground 
states |0, 0), Itt, 0), |0,7r), |7r, tt) are linearly independent. 
To obtain the entanglement spectrum for a given ground 
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FIG. 2: Above: Renyi entropy S2 for the state 1$) — 
coa{(j))\Q,Qi) -f sm(0) |0,7r) for different values of the gap 
A that enters the CSL wavefunction. The dotted curve 
is the theoretical prediction in the thermodynamic limit: 
S2 = C(A) - ln(cos''(0) + sin^{ct))) where we C(A) is a non- 
universal constant (here we chose C = 1.35 so as to fa- 
cilitate the comparison with the numerical data). Further, 
S2,max — S2,min ~ 0.67 independent of A which is close to the 
expected result of ln(2) in the thermodynamic limit. Below: 
The difference between 52 for two sets of two different A's as 
a function of the angle 4>. The difference 5*2 (Ai) — 5*2 (A2), is 
nearly (j) independent in agreement with the theoretical pre- 
diction, since in it equals S2,iocai{^i)- S2,iocai{^2) (see Eq|6j. 



state, we divide the system into two halves A and B, each 
of size 3x3 and express the corresponding wavefunction 
as 



\a,) \bj) 



(7) 



where \ai) and \bi) are the basis states for the region 
A and B. If the regions A and B have A^/2 spins each, 
then is the ij'th element of a 2^/^ x 2^/^ matrix. 
The Schmidt decomposition of the wavefunction |^) (see 
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0,0) 



FIG. 3: Eigenvalues Aq of the density matrix corresponding 
to the two states |0,0) and |7r,0). Since very large number 
of eigenvalues have a very small value, only few of them are 
visible. The spectrum for the two states shown is identical 
within computer precision. 



Eql2]) is obtained by performing the Singular Value De- 
composition (SVD) of the matrix ipij numerically. 

For the system size considered (6x3 lattice) , the MESs 
are found to correspond to the states |Si) = |0, 0) and 
IS2) = |0,7r) (Figl2]) . Because of finite size effects, the 
orthogonality of these states is not perfect for this sys- 
tem size with (0,7r|0,0) w 0.03 - 0.07 for the range of 
mean-field gap A considered. Despite this, we find that 
the entanglement spectrum of these two states is identi- 
cal (Fig|3]) in agreement with EqlS] We do not resolve the 
momentum along the cut for the Schmidt states, since it 
is not needed for our purposes. Further, we find that the 



Schmidt states 



corresponding to same value of 



the density-matrix eigenvalue Aq, are orthogonal in agree- 
ment with EqlH However, because of lack of the com- 
plete orthogonality of the two MES, the Schmidt states 
for non-identical eigenvalues arc found to be not always 
orthogonal. 

FigE] shows the Renyi entropy S'2 for the linear com- 
bination 1$) = cos(0) |0,0) + sin{<j))\0,Tr) for different 
values of the mean-field d-wave BCS gap A, that enters 
as a parameter in the CSL wavefunction. We notice that 
the difference S2,iocaii^i, 4>) - S2,iocai{^2, (t>) is nearly in- 
dependent of (j) in agreement with Eq|6l the difference 
depends only on Siocai{^i) - Siocai{^2), thereby being 
only a function of Ai,A2. Perhaps most interestingly, 
Eq|6] predicts that for the CSL state, S2,max — S2,min = 
'^Itopo = ln(2) « 0.69. We find that the numerical value 
of the same quantity equals ~ 0.67, in close agreement 
with the theory, and remains essentially constant as A is 
changed. We note that though the curve S2{4') in Figl3] 
docs not agree with the theoretical prediction in the ther- 



modynamic limit completely because of the aforemen- 
tioned finite size effects , the overall agreement it is still 
remarkable with a system size of only 18 sites. 

Since the matrix ij} is non-sparse, the numerical cost of 
performing SVD on a 2^^"^ x 2^/^ is nearly equal to the 
cost of exact diagonalizing the whole system for a local 
Hamiltonian. Therefore, our scheme can be implemented 
exactly on system sizes upto iV = 30 — 40 sites for spin 
S = 1/2. Of course, the detection of topological order 
in such sizes will be contingent on whether the correla- 
tion length of various local operators ^ <^ L, the linear 
length of the system. As mentioned in the introduction, 
it is rather fortuitous that one of the best current meth- 
ods of choice to solve strongly correlated systems, namely 
DMRG, happens to provide one with the entanglement 
spectrum of the ground state wavcfunction(s). Given the 
hints of topological order in KAH with two nearly degen- 
erate ground states and the short correlation length of 
one- two unit cells observed in Ref.jl^, it may be worth- 
while to apply the method described in this paper to 
KAH using DMRG and exact diagonalization as well. 

Finally, some the considerations developed above 
are also applicable to the so-called "quantum-ordered" 
phases, namely, phases with emergent gaplcss gauge 
bosons that may be coupled to bosonic/fcrmionic matter 
0. As an example, consider the U(l) spin-liquid in three 
dimensions (see, e.g. Ref.[38|). The low-energy spectrum 
of this spin-liquid contains analog of topological degener- 
ate states whose gap scale as ~ 1/L with the system size 
L. One might think that this renders them difficult to 
distinguish from the excitations corresponding to fluctu- 
ations of emergent photons whose gap also scales as 1/L. 
However, based on the discussion in this paper, entan- 
glement properties can, in fact, distinguish states that 
correspond to these two distinct types of excitations. In 
particular, for a given entanglement cut, using the states 
that arc analog of topological degeneracy, one may again 
construct a set of minimum entropy states that will sat- 
isfy Eqini while such a relation will generically not hold 
for the other low-lying states. 

To summarize, we described a method that can be 
used to detect topological order from the Schmidt de- 
composed degenerate ground states. The essential part 
consists in noticing that the non-universal leading term 
in the entanglement entropy is identical for all degener- 
ate ground states and thus can be canceled out so as to 
extract the subleading universal constant, whose form is 
dictated by EqlBJ when the boundary of the bipartition 
is non-contractible. We also demonstrated the method 
by performing exact Schmidt decomposition of a topo- 
logically ordered state, the chiral spin liquid state. We 
believe that this method can be put to practical use in 
numerical methods such as DMRG and exact digonaliza- 
tion. 
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